Cell sorting by active forces in a phase-field model of cell monolayers.

Cell sorting, the segregation of cells with different properties into distinct domains, is a key phenomenon in biological processes such as embryogenesis. We use a phase-field model of a confluent cell layer to study the role of activity in cell sorting. We find that a mixture of cells with extensile or contractile dipolar activity, and which are identical apart from their activity, quickly sort into small, elongated patches which then grow slowly in time. We interpret the sorting as driven by the different diffusivity of the extensile and contractile cells, mirroring the ordering of Brownian particles connected to different hot and cold thermostats. We check that the free energy is not changed by either partial or complete sorting, thus confirming that activity can be responsible for the ordering even in the absence of thermodynamic mechanisms.

Shape-Tension Coupling Produces Nematic Order in an Epithelium Vertex Model.

We study the vertex model for epithelial tissue mechanics extended to include coupling between the cell shapes and tensions in cell-cell junctions. This coupling represents an active force which drives the system out of equilibrium and leads to the formation of nematic order interspersed with prominent, long-lived +1 defects. The defects in the nematic ordering are coupled to the shape of the cell tiling, affecting cell areas and coordinations. This intricate interplay between cell shape, size, and coordination provides a possible mechanism by which tissues could spontaneously develop long-range polarity through local mechanical forces without resorting to long-range chemical patterning.

Active nematics with deformable particles.

The hydrodynamic theory of active nematics has been often used to describe the spatio-temporal dynamics of cell flows and motile topological defects within soft confluent tissues. Those theories, however, often rely on the assumption that tissues consist of cells with a fixed, anisotropic shape and do not resolve dynamical cell shape changes due to flow gradients. In this paper we extend the continuum theory of active nematics to include cell shape deformability. We find that circular cells in tissues must generate sufficient active stress to overcome an elastic barrier to deforming their shape in order to drive tissue-scale flows. Above this threshold the systems enter a dynamical steady-state with regions of elongated cells and strong flows coexisting with quiescent regions of isotropic cells.

Geometrical control of interface patterning underlies active matter invasion.

Interaction between active materials and the boundaries of geometrical confinement is key to many emergent phenomena in active systems. For living active matter consisting of animal cells or motile bacteria, the confinement boundary is often a deformable interface, and it has been unclear how activity-induced interface dynamics might lead to morphogenesis and pattern formation. Here, we studied the evolution of bacterial active matter confined by a deformable boundary. We found that an ordered morphological pattern emerged at the interface characterized by periodically spaced interfacial protrusions; behind the interfacial protrusions, bacterial swimmers self-organized into multicellular clusters displaying +1/2 nematic defects. Subsequently, a hierarchical sequence of transitions from interfacial protrusions to creeping branches allowed the bacterial active drop to rapidly invade surrounding space with a striking self-similar branch pattern. We found that this interface patterning is geometrically controlled by the local curvature of the interface, a phenomenon we denote as collective curvature sensing. Using a continuum active model, we revealed that the collective curvature sensing arises from enhanced active stresses near high-curvature regions, with the active length scale setting the characteristic distance between the interfacial protrusions. Our findings reveal a protrusion-to-branch transition as a unique mode of active matter invasion and suggest a strategy to engineer pattern formation of active materials.

Phase Separation Driven by Active Flows.

We extend the continuum theories of active nematohydrodynamics to model a two-fluid mixture with separate velocity fields for each fluid component, coupled through a viscous drag. The model is used to study an active nematic fluid mixed with an isotropic fluid. We find microphase separation, and argue that this results from an interplay between active anchoring and active flows driven by concentration gradients. The results may be relevant to cell sorting and the formation of lipid rafts in cell membranes.

Collective rotational motion of freely expanding T84 epithelial cell colonies.

Coordinated rotational motion is an intriguing, yet still elusive mode of collective cell migration, which is relevant in pathological and morphogenetic processes. Most of the studies on this topic have been carried out on epithelial cells plated on micropatterned substrates, where cell motion is confined in regions of well-defined shapes coated with extracellular matrix adhesive proteins. The driver of collective rotation in such conditions has not been clearly elucidated, although it has been speculated that spatial confinement can play an essential role in triggering cell rotation. Here, we study the growth of epithelial cell colonies freely expanding (i.e. with no physical constraints) on the surface of cell culture plates and focus on collective cell rotation in such conditions, a case which has received scarce attention in the literature. One of the main findings of our work is that coordinated cell rotation spontaneously occurs in cell clusters in the free growth regime, thus implying that cell confinement is not necessary to elicit collective rotation as previously suggested. The extent of collective rotation was size and shape dependent: a highly coordinated disc-like rotation was found in small cell clusters with a round shape, while collective rotation was suppressed in large irregular cell clusters generated by merging of different clusters in the course of their growth. The angular motion was persistent in the same direction, although clockwise and anticlockwise rotations were equally likely to occur among different cell clusters. Radial cell velocity was quite low as compared to the angular velocity, in agreement with the free expansion regime where cluster growth is essentially governed by cell proliferation. A clear difference in morphology was observed between cells at the periphery and the ones in the core of the clusters, the former being more elongated and spread out as compared to the latter. Overall, our results, to our knowledge, provide the first quantitative and systematic evidence that coordinated cell rotation does not require a spatial confinement and occurs spontaneously in freely expanding epithelial cell colonies, possibly as a mechanism for the system.

Active Forces in Confluent Cell Monolayers.

We use a computational phase-field model together with analytical analysis to study how intercellular active forces can mediate individual cell morphology and collective motion in a confluent cell monolayer. We explore the regime where intercellular forces dominate the tissue dynamics, and polar forces are negligible. Contractile intercellular interactions lead to cell elongation, nematic ordering, and active turbulence characterized by motile topological defects. Extensile interactions result in frustration, and perpendicular cell orientations become more prevalent. Furthermore, we show that contractile behavior can change to extensile behavior if anisotropic fluctuations in cell shape are considered.

Steering self-organisation through confinement.

Self-organisation is the spontaneous emergence of spatio-temporal structures and patterns from the interaction of smaller individual units. Examples are found across many scales in very different systems and scientific disciplines, from physics, materials science and robotics to biology, geophysics and astronomy. Recent research has highlighted how self-organisation can be both mediated and controlled by confinement. Confinement is an action over a system that limits its units' translational and rotational degrees of freedom, thus also influencing the system's phase space probability density; it can function as either a catalyst or inhibitor of self-organisation. Confinement can then become a means to actively steer the emergence or suppression of collective phenomena in space and time. Here, to provide a common framework and perspective for future research, we examine the role of confinement in the self-organisation of soft-matter systems and identify overarching scientific challenges that need to be addressed to harness its full scientific and technological potential in soft matter and related fields. By drawing analogies with other disciplines, this framework will accelerate a common deeper understanding of self-organisation and trigger the development of innovative strategies to steer it using confinement, with impact on, e.g., the design of smarter materials, tissue engineering for biomedicine and in guiding active matter.

Activity-driven tissue alignment in proliferating spheroids.

We extend the continuum theory of active nematic fluids to study cell flows and tissue dynamics inside multicellular spheroids, spherical, self-assembled aggregates of cells that are widely used as model systems to study tumour dynamics. Cells near the surface of spheroids have better access to nutrients and therefore proliferate more rapidly than those in the resource-depleted core. Using both analytical arguments and three-dimensional simulations, we find that the proliferation gradients result in flows and in gradients of activity both of which can align the orientation axis of cells inside the aggregates. Depending on environmental conditions and the intrinsic tissue properties, we identify three distinct alignment regimes: spheroids in which all the cells align either radially or tangentially to the surface throughout the aggregate and spheroids with angular cell orientation close to the surface and radial alignment in the core. The continuum description of tissue dynamics inside spheroids not only allows us to infer dynamic cell parameters from experimentally measured cell alignment profiles, but more generally motivates novel mechanisms for controlling the alignment of cells within aggregates which has been shown to influence the mechanical properties and invasive capabilities of tumors.

Viscoelastic confinement induces periodic flow reversals in active nematics.

We use linear stability analysis and hybrid lattice Boltzmann simulations to study the dynamical behavior of an active nematic confined in a channel made of viscoelastic material. We find that the quiescent, ordered active nematic is unstable above a critical activity. The transition is to a steady flow state for high elasticity of the channel surroundings. However, below a threshold elastic modulus, the system produces spontaneous oscillations with periodic flow reversals. We provide a phase diagram that highlights the region where time-periodic oscillations are observed and explain how they are produced by the interplay of activity and viscoelasticity. Our results suggest experiments to study the role of viscoelastic confinement in the spatiotemporal organization and control of active matter.

Activity gradients in two- and three-dimensional active nematics.

We numerically investigate how spatial variations of extensile or contractile active stress affect bulk active nematic systems in two and three dimensions. In the absence of defects, activity gradients drive flows which re-orient the nematic director field and thus act as an effective anchoring force. At high activity, defects are created and the system transitions into active turbulence, a chaotic flow state characterized by strong vorticity. We find that in two-dimensional (2D) systems active torques robustly align +1/2 defects parallel to activity gradients, with defect heads pointing towards contractile regions. In three-dimensional (3D) active nematics disclination lines preferentially lie in the plane perpendicular to activity gradients due to active torques acting on line segments. The average orientation of the defect structures in the plane perpendicular to the line tangent depends on the defect type, where wedge-like +1/2 defects align parallel to activity gradients, while twist defects are aligned anti-parallel. Understanding the response of active nematic fluids to activity gradients is an important step towards applying physical theories to biology, where spatial variations of active stress impact morphogenetic processes in developing embryos and affect flows and deformations in growing cell aggregates, such as tumours.

Active Extensile Stress Promotes 3D Director Orientations and Flows.

We use numerical simulations and linear stability analysis to study an active nematic layer where the director is allowed to point out of the plane. Our results highlight the difference between extensile and contractile systems. Contractile stress suppresses the flows perpendicular to the layer and favors in-plane orientations of the director. By contrast extensile stress promotes instabilities that can turn the director out of the plane, leaving behind a population of distinct, in-plane regions that continually elongate and divide. This supports extensile forces as a mechanism for the initial stages of layer formation in living systems, and we show that a planar drop with extensile (contractile) activity grows into three dimensions (remains in two dimensions). The results also explain the propensity of disclination lines in three dimensional active nematics to be of twist type in extensile or wedge type in contractile materials.

Flow transitions and length scales of a channel-confined active nematic.

We perform lattice Boltzmann simulations of an active nematic fluid confined in a two-dimensional channel to study the range of flow states that are stabilised by the confinement: unidirectional flow, oscillatory flow, the dancing state, localised active turbulence and fully-developed active turbulence. We analyse the flows in Fourier space, and measure a range of different length scales which describe the flows. We argue that the different states occur as a result of flow instabilities inherent to the system. As a consequence the characteristic length scale for oscillatory flow, the dancing state and localised active turbulence is set by the channel width. Fully-developed active turbulence occurs only when the channel width is larger than the intrinsic, active length scale of the bulk fluid. The results clarify why the activity number is a control parameter for the flow transitions.

Coupling Turing stripes to active flows.

We numerically solve the active nematohydrodynamic equations of motion, coupled to a Turing reaction-diffusion model, to study the effect of active nematic flow on the stripe patterns resulting from a Turing instability. If the activity is uniform across the system, the Turing patterns dissociate when the flux from active advection balances that from the reaction-diffusion process. If the activity is coupled to the concentration of Turing morphogens, and neighbouring stripes have equal and opposite activity, the system self organises into a pattern of shearing flows, with stripes tending to fracture and slip sideways to join their neighbours. We discuss the role of active instabilities in controlling the crossover between these limits. Our results are of relevance to mechanochemical coupling in biological systems.

Submersed micropatterned structures control active nematic flow, topology, and concentration.

Coupling between flows and material properties imbues rheological matter with its wide-ranging applicability, hence the excitement for harnessing the rheology of active fluids for which internal structure and continuous energy injection lead to spontaneous flows and complex, out-of-equilibrium dynamics. We propose and demonstrate a convenient, highly tunable method for controlling flow, topology, and composition within active films. Our approach establishes rheological coupling via the indirect presence of fully submersed micropatterned structures within a thin, underlying oil layer. Simulations reveal that micropatterned structures produce effective virtual boundaries within the superjacent active nematic film due to differences in viscous dissipation as a function of depth. This accessible method of applying position-dependent, effective dissipation to the active films presents a nonintrusive pathway for engineering active microfluidic systems.

Activity pulses induce spontaneous flow reversals in viscoelastic environments.

Complex interactions between cellular systems and their surrounding extracellular matrices are emerging as important mechanical regulators of cell functions, such as proliferation, motility and cell death, and such cellular systems are often characterized by pulsating actomyosin activities. Here, using an active gel model, we numerically explore spontaneous flow generation by activity pulses in the presence of a viscoelastic medium. The results show that cross-talk between the activity-induced deformations of the viscoelastic surroundings and the time-dependent response of the active medium to these deformations can lead to the reversal of spontaneously generated active flows. We explain the mechanism behind this phenomenon based on the interaction between the active flow and the viscoelastic medium. We show the importance of relaxation time scales of both the polymers and the active particles and provide a phase space over which such spontaneous flow reversals can be observed. Our results suggest new experiments investigating the role of controlled pulses of activity in living systems ensnared in complex mircoenvironments.

Investigating the nature of active forces in tissues reveals how contractile cells can form extensile monolayers.

Actomyosin machinery endows cells with contractility at a single-cell level. However, within a monolayer, cells can be contractile or extensile based on the direction of pushing or pulling forces exerted by their neighbours or on the substrate. It has been shown that a monolayer of fibroblasts behaves as a contractile system while epithelial or neural progentior monolayers behave as an extensile system. Through a combination of cell culture experiments and in silico modelling, we reveal the mechanism behind this switch in extensile to contractile as the weakening of intercellular contacts. This switch promotes the build-up of tension at the cell-substrate interface through an increase in actin stress fibres and traction forces. This is accompanied by mechanotransductive changes in vinculin and YAP activation. We further show that contractile and extensile differences in cell activity sort cells in mixtures, uncovering a generic mechanism for pattern formation during cell competition, and morphogenesis.

Memory effects, arches and polar defect ordering at the cross-over from wet to dry active nematics.

We use analytic arguments and numerical solutions of the continuum, active nematohydrodynamic equations to study how friction alters the behaviour of active nematics. Concentrating on the case where there is nematic ordering in the passive limit, we show that, as the friction is increased, memory effects become more prominent and +1/2 topological defects leave increasingly persistent trails in the director field as they pass. The trails are preferential sites for defect formation and they tend to impose polar order on any new +1/2 defects. In the absence of noise and for high friction, it becomes very difficult to create defects, but trails formed by any defects present at the beginning of the simulations persist and organise into parallel arch-like patterns in the director field. We show aligned arches of equal width are approximate steady state solutions of the equations of motion which co-exist with the nematic state. We compare our results to other models in the literature, in particular dry systems with no hydrodynamics, where trails, arches and polar defect ordering have also been observed.